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Abstract 

Motivated by conventional gauge theories, we consider a theory of 
gravity in which the Einstein-Hilbert action is replaced by a term that is 
quadratic in the Riemann tensor. We focus on cosmological solutions to 
the field equations in flat, open and closed universes. The gravitational 
' action is scale invariant, so the only matter source considered is radiation. 

, The theory can also accommodate isotropic torsion and this generically 

£S| ■ removes singularities from the evolution equations. For general initial 

conditions the Hubble parameter H(i) is driven in a seemingly chaotic 
fashion by torsion to produce irregularly occuring inflationary regions. In 
■ the absence of torsion, the theory reproduces the standard cosmological 

solutions of a simple big bang model. A satisfying feature is that a cos- 
mological constant arises naturally as a constant of integration, and does 
not have to be put into the Lagrangian by hand. 



in 



Or 



(N 



> , 

^ ■ 1 Introduction 

O 

Q\ • Many authors have argued that gravity can be formulated as a gauge the- 

ory ^ |2 Ej . But there is a fundamental difference between the Lagrangians 
of the electroweak and strong interactions, and the Einstein-Hilbert action of 
gravity. The former are all quadratic in the field strength tensor, whereas the 
latter is only linear. It is therefore of interest to consider the consequences of a 



theory of gravity based on a Lagrangian that is quadratic in the Riemann tensor. 



1 ■ This is the theory we consider here. We retain the gauge-theoretic framework 

of the Einstein-Cartan-Kibble-Sciama theory 0], so the tetrad (vierbein) and 
connection are varied independently. Furthermore, no derivatives of the vier- 
bein enter the Lagrangian, so that its variation is straightforward and gives the 
j_j ■ quadratic analog of the Einstein equations. 

A problematic feature of the gravitational Lagrangian we consider is that it 
is scale invariant. This implies that the only matter stress-energy tensors that 
can couple to gravity are trace-free. Ultimately this symmetry would have to 
be broken in any physically realistic theory. But here we circumvent the prob- 
lem by only considering cosmological solutions where the matter is described 
by radiation (P = p/3). A more attractive property of a quadratic Riemann 
Lagrangian is that torsion can be generated without an accompanying source of 
spin |S] . This is explicitly prohibited in Ricci scalar Lagrangians JU ■ The torsion 
is not restricted in type, and can take the form of an isotropic field appropriate 
for cosmology. (For a review of cosmological models with anisotropic torsion, 
see 0].) 
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Here we construct the full cosmological field equations in a straightforward 
manner, allowing for the presence of a torsion field. The symmetries in the 
Lagrangian allow the equations to be expressed in a variety of ways, and we 
show that one set of gauge choices leads to a simple complex-valued differential 
equation. This equations can be solved explicitly in terms of elliptic functions, 
though the accompanying evolution looks almost chaotic. We find in particular 
that the torsion couples to the Hubble parameter to produce inflationary regimes 
in a quasi-periodic fashion. While the present day experimental upper bound on 
background cosmological torsion is very small, on the order of 10 -30 GeV jEllZh 
our model shows that it could have been very large during some earlier epoch, 
and could become large again in the future. 

A further attractive feature of the theory presented here is that many of the 
standard cosmological models (coupled to radiation) are contained as special 
cases. With zero torsion the equations are solved by a standard radiation- 
dominated cosmology. This set-up illustrates neatly how the inclusion of a 
small amount of torsion removes the big bang singularity by translating it into 
the complex plane. Minkevich |S] and Faria-Busto j) have also argued that a 
quadratic theory can remove the big bang singularity. Furthermore, we are also 
able to find solutions that have a cosmological constant present, despite the fact 
that is has not been put in by hand in the original Lagrangian. This happens 
because the equations are now of higher order, and a cosmological constant 
arises naturally as a constant of integration. 

The paper begins with a derivation of the field equations. Much of the 
complexity of this derivation is avoided by working at the level of the action. 
We then study some of the general properties of solutions, before considering 
some special cases. Unless otherwise stated we work in natural units where 
G = c= 1. 



2 The field equations 

Our starting point is the action integral 

S = J d 4 x^g \^R a ^»Rj fil/ + kL^ (1) 

where R a is the Riemann tensor, C m is the matter Lagrangian, and k = 8n. 
The gravitational degrees of freedom in this theory are a tetrad h^ a and the 
connection ui fl ab . Throughout we use Greek letters for coordinate components, 
and Latin for tetrad components. The Riemann tensor is constructed from the 
connection by 

nab o ab a ab , , , ac b ac b fry\ 

and this can be thought of as the field strength of the connection field. 

Though the theory discussed here is clearly different from general relativity, 
it is still often convenient to talk in terms of Riemannian geometry and torsion. 
The metric is defined in the expected manner by 

= Vabh^hJ 3 '. (3) 
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The metric defines a Christoffcl connection , and from this we can write 

= h a a (T x „ + K x ^)h f \ - d^h [3 a h f \. (4) 

This equation defines the contorsion tensor K x ^ v in terms of the pieces of uj^ a b 
that would not be present in pure, torsion-free, general relativity. 

In this paper we are interested in cosmological solutions. As the underlying 
metric still gives rise to many of the observational consequences of the theory, the 
standard cosmological arguments apply and we can restrict to a homogeneous, 
isotropic line element of the form 



(ds) 2 = dt 2 - a{tf 



dr 2 



1 — kr 2 



r 2 {d0 2 + sin 2 9 d<p 2 ) 



(5) 



Here k € {±1, 0} sets the curvature, a(t) is the scale factor, and the bar on <p is 
present as the (p symbol will be used elsewhere. A suitable tetrad for this metric 
is given by 

(I \ 

a(l-kr 2 )-i 
ar 

\0 arsinOj 

The most general spin connection uj^ a b consistent with isotropy and homo- 
geneity has components 



h a 



(6) 



Wioi = 
W21 2 = yl - kr 2 , 



VT— kr 2 



^123 — — 



^312 = — rip sin( 



^231 
^313 



Vl — kr 2 ' 
-rip, 

yl - kr 2 sin ( 



^202 = r<p, 

^303 = r0sin6», 
W323 = cos0, (7) 



together with further terms generated by the antisymmetry in the last two in- 
dices. Here <p and ip are scalar functions of time. The theory therefore contains 
three gravitational degrees of freedom — a(t), <p(t) and ip(t). It is now a straight- 
forward matter to compute the Riemann tensor and, in terms of a tetrad basis, 
we find the following nonzero components: 



R°joi = 


-Sji 
a 


lx - ijk 


R l 0k0 = 


— d k 


R % 0jk 


R l j0k — 


ip 

^ijk 

a 


R l jkl 



20V 
a 2 

20V 



2 £i 3 k 
-lP 2 +k) 

~a~ 2 



a 

k2 _ „/,2 



(S l k Sji - SjSjk) 



(8) 



(together with further terms generated by the antisymmetry in the last two 
indices). 

The spatial integrals are now irrelevant to the full action integral, and the 
cosmological equations are generated by the reduced action 

S = j dt (a(tf - V 2 ) + ^(0 4 - 6<p 2 ip 2 + V 4 + 2fc0 2 - 2fcV 2 + k 2 ) + ^a 3 p) , 

(9) 
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where we have assumed that the matter content is pure radiation (p = p/3). 
This implies that the matter content does not contain any direct coupling to 
the spin connection. 

Scale invariance is still present in this action, though in the slightly modified 
form of a combined rescaling and reparameterisation. The effect of these is that 
the scale factor a transforms as the inverse of an einbein, 

a(t) » a'(t) = J^a(i'), (10) 

with and ip transforming as scalar fields, cf>(t) i— > 4>'(t) — 4>(t'). The result of 
this symmetry is that the scale factor a is essentially an arbitrary function and, 
after variation, can be chosen in any convenient manner. Two choices of a are 
particularly attractive. The first is the 'cosmological' gauge choice, for which 

<t> = a, (11) 

where the overdot represents differentiation with respect to t. This is the gauge 
that most closely resembles a standard FRW cosmology. The evolution param- 
eter t measures cosmic time, and the equations can be expressed in terms of the 
Hubble function H(t), where 

H=-- (12) 
a 

The second gauge of interest is the 'conformal' gauge in which we set a = 1. In 
this case the evolution parameter t measures conformal time, which we denote 
by r\. The conformal gauge turns out to be most convenient for solving the field 
equations. 

3 Equations of motion and their general solu- 
tion 

The Euler-Lagrange equations for <p and ip from the action integral J^J) are 

dt(a^) = -U 3 -3^ 2 +k(f)) (13) 
a 

-d t (aip) = -(-30V + V> 3 - kip) (14) 
a 

The a variation produces a form of Hamiltonian constraint equation, 

a(0 2 - ifj 2 ) ^-((f)*- 60 2 t/> 2 + 4> A + 2k(j) 2 - 2ki> 2 + k 2 ) - na 3 p. (15) 
a 

An immediate consequence of these equations is that we recover the expected 
conservation law for radiation, 

a 4 p = constant. (16) 

In the absence of torsion %p = and, working in the cosmological gauge, the 
dynamics reduces to the pair of equations 

k 

H + AHH = 2H— (17) 
a 1 

(H + H 2 ) 2 = (h 2 + - Kp . (18) 
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We return to these equations in section 0] 

The full equations of motions are simplest to analyse in the conformal gauge. 
Setting a = 1, so that all derivatives are with respect to conformal time, we 
arrive at the pair of equations 

(/}" = 2<j) 3 - 6cj)ip 2 + 2k(t) (19) 
ip" = -60 2 i/> + 2^ 3 - 2kifj. (20) 

Henceforth we use dashes to denote derivatives with respect to conformal time, 
and retain overdots for derivatives with respect to cosmic time. We now intro- 
duce the complex variable 

Z = </> + iip (21) 

so that the equations of motion are contained in the single complex-valued 
differential equation 

Z" = 2Z 3 + 2kZ. (22) 
This equation has a trivial first integral, 

Z' 2 = Z 4 + 2kZ 2 + A (23) 

where Aq is a complex- valued constant. The real part of Aq is determined by 
the conserved energy density a 4 p and the curvature via 

Rc(^l ) = k 2 - na 4 p. (24) 

The imaginary part of Aq can similarly be thought of a source term for the 
torsion. The first integral (|23[l illustrates that the system is essentially Hamil- 
tonian, but the complex structure drives a mixed signature in momentum space 
(going as tfi' 2 — ip' 2 ), which has some unusual consequences. 

From the first integral form of the equations of motion, the full solution can 
be written down as 

r dz 

V ~J (Z 4 + 2fcZ 2 + A )i/2- (25) 

This integral can be solved analytically in terms of the Weierstrass elliptic func- 
tion p(rf) = p(rj; 52,33) |10l 1111 IT2*] . This function requires two constants, 32 
and 33. These are defined in terms of the quartic polynomial 

f(Z) = Z 4 + 2kZ 2 + A . (26) 

From this we compute the pair of elliptic invariants (72 and g%, obtaining 

g 2 =A + 12k 2 (27) 
.93 = 4fcA . (28) 

The integral 

r z dZ' 

now has the general solution 



z = z + yggoVfa) + into) (pM ~ jAgo)) + jif(Z )f"(Z ) 
2(p[rj) - ±f"(Z )) 2 ±f(Z )fW(Z Q ) 
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Figure 1 : Evolution of Z without radiation. The plot shows the real (solid) and 
imaginary (dashed) components of Z{rj) as a function of conformal time r\. For 
these plots we have f(Z) = Z 4 — i, Zq = 1 — i and the sign for the initial value 
of Z' is chosen so that is decreasing. Also highlighted are the regions in which 
inflation is occurring, where < <f)' < <fi 2 . The magnitude of the inflation is 
desribed by <f>' . 



The Weierstrass p(?7)-function has a second order pole at n = 0, and the deriva- 
tive of the Weierstrass function, p'{rj), is itself an elliptic function with pole of 
order 3 at rj = 0. The constant Zq is arbitrary, and it is sometimes convenient 
to chose it at one of the zeros of f(Z). With f(Z ) — the solution simplifies 
to 

Z = Z Q + \f{Z ) L(r,; g2, gS) - ^/"(^o)) • (31) 

Given generic initial conditions, the behavior of the general solution (|30|l is 
extemely complicated, as the following plots illustrate. Figures H and [2 show 
sample behaviour both without and then with matter. We can generally charac- 
terize the general solution as describing periods of fairly stable behaviour, which 
are quasi-periodically driven by the torsion to produce very large, but generally 
finite, peaks. 

The plots Hand [3 also highlight the regions where inflation is occurring. To 
define the regions of inflation we must first recall how to switch between the 
conformal and cosmological gauges. Given (f>(rj) we can introduce a scale factor 
via equation lfTT|) , the solution of which is 

a(*7)=expQ^(y)*/)- (32) 
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Figure 2: Evolution of Z with radiation. A similar graph to figure ^ but now 
with matter present. This in incorporated by setting Ao = 1 — i. We keep the 
initial value of Zq = 1 — i. 



With a computed, we then compute the Hubble function by 

If necessary, H (rj) can be converted back to a function of cosmic time, though 
it is easier to compute H(t) from the differential equations in the cosmic gauge. 
Typical behaviour for H{t) is shown in figure 

In the theory as it stands it makes no sense to distinguish one gauge as 
being preferred over another, as we have no way of measuring H, for example. 
One would have to find some way of breaking the conformal invariance and 
including matter with mass. However, it is still interesting to compute whether 
the solutions could represent an inflating universe. The standard condition 
for inflation is d > 0, and we also impose the further restriction that H < 
0. Combined, these produce the restriction < — H < H 2 . In terms of the 
conformal gauge, the equivalent condition is 

< <j)' < <j) 2 . (34) 

The regions for which this occurs are marked on figures and [5] 

A further way of viewing the dynamics is shown in figure 0] This plot 
shows the trajectory of Z(t) in complex space over a finite time interval. The 
plot further illustrates how the complex dynamical evolution can arise from a 
simple complex ordinary differential equation. The plot concentrates on the fine 
structure around the origin and does not show the large excursions that occur 
aperiodically. 
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Figure 3: Hubble function H as a function of cosmic time t. The initial condi- 
tions are chosen to match those of figure JIJ. Also plotted is the torsion Q(t), 
where Q = ip/a. Inflationary regions are also highlighted. 




Figure 4: Complex space trajectory of Z . The plot shows the trajectory in 
complex space of Z[rj) = <j> + itp over a finite interval of conformal time. The 
plot is for f(Z) = Z 4 — i with initial condition Zq = 1 — i. The plot illustrates 
how the trajectory of Z(r\) generally stays in the vicinity of the origin, but 
occasionally moves far away. The horiztonal line illustrates one such excursion. 
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4 Special cases 



In this section we consider some of the special case solutions to the general 
cosmological equations. First, returning to the torsion-free equation (|17|l and 
integrating we find that 

H + 2H 2 = -4+ co (35) 
cr 

where cq is a constant of integration. It then follows from equation (|18|l that 
the density is given by 

^P = H 2 + ^ + ^. (36) 
Acq 2 a z 

But, following a redefinition of p, this is precisely the standard FRW equation 
with a cosmological constant Co/2. Indeed, in order to include matter into this 
model, a term representing the cosmological constant is forcedl This is clear 
from the fact that cq = implies p = 0. 

If we continue to specialise by further setting k and p to zero, the cosmolog- 
ical equations reduce to 

H + H 2 = ±H 2 . (37) 

This has two solutions. The positive root corresponds to constant H, repre- 
senting a de Sitter space. The negative root has H = —2H 2 and is solved 
by H = \/(2t). This corresponds to a simple big bang model. So we see that 
many standard cosmological solutions are contained in the quadratic Lagrangian 
framework. 

It is instructive to see how we recover the simple big bang solution H = 
l/(2i) from the general solution In the conformal gauge, this solution is 

characterised by 

z =h ^ 

To recover this solution we first set f(Z) = Z , which removes the matter and 
curvature contributions. The general solution now reduces to 1)30(1 to 

Zjp'jfj) + 2Zlp(r 1 ) 
Z - Z ° + 2p*(„) - 2Z 2 P (V) • ^ 

(Selecting the opposite sign for the square root in this expression only changes 
the overall sign in the result.) With 52 = S3 = we have 

p{v) = r\ (40) 

and it follows that 

7 -^—T- (41) 



V + z o 

By choosing Zq 1— > 00 we recover the big bang model. An interesting effect 
occurs if we include some initial torsion by setting by setting Zq = —i/b. Now 
the general solution becomes 

„ ill — ib , . » 

z = 7 9 42 
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Figure 5: Repeated big bang and big crunch scenario without torsion. Initial 
values are Aq = — 1 and Zq = 1, corresponding to elliptic invariants of 32 = — 1 
and 33 = 0. This is known as the pseudolcmniscate case. 



We see now that the initial big bang singularity has been removed by the torsion 
and the evolution is entirely finite. In effect, the torsion term has moved the 
singularity to a complex value of the conformal time. 

We can find further interesting models within the general solution. Choosing 
f(Z) = Z 4 + 1 and Zq = 1 then the solution produces big bang and big crunch 
singularities in a periodic manner (see figure EJ). For this configuration, 32 = 1 
and <7i = 0. This corresponds to one of three special cases of the elliptic invari- 
ants that has been given a special name, this one being the pseudolemniscate 
case The pseudolemniscate case has half-periods L(±l + i)/4 where L is 
the lemniscate constant L/4 = 1.311028 . . .. 



5 Conclusions 

We have shown that for a Lagrangian that is quadratic in the Riemann tensor 
the cosmological equations can be solved analytically in terms of the elliptic 
Weierstrass function p(j]). The model can only couple to scale-invariant fields, 
and the only matter considered is isotropic radiation. Furthermore, we have 
shown that these solutions reduce to standard cosmological solutions in various 
limits. This model is particularly interesting in that it produces many periods 
of inflation for a wide range of initial conditions. But in order for such a model 
to be physical it is necessary to break conformal invariance and introduce mass, 
which will form the subject of future research. 
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